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ABSTRACT 


In the present work, we write a brief exposition of the Jefimenko’s theory of gravitation. This theory arose 
from the analogy between the laws of gravitation and electromagnetism, that is, exist a second gravitational 
field called cogravitational field, analogous to the magnetic field. We introduce a new system of units called 
Gravitational Gaussian System (GGS). This system allows us write the equations of gravitation in a simple 
form to solve them. Using the Jefimenko’s equations of gravitation, we obtain the wave equations for 
gravitational and cogravitational fields and we find wave solutions. We demonstrate there are 


configurations of the gravitodynamical field (that is, the set of gravitational and cogravitacional fields) in 
form of co gravitational field spheres and gravitational field rings. This phenomenon must be an analogue of 
the ball lightning in the electromagnetic field, but in this case the cogravitational field spheres serves as 
containers of matter (it could be a gas). We analyze how this configuration acts on the particles inside the 
spheres, and we investigate the physical properties of such configurations, namely, how behaves the density 


of energy and the Pointing vector of this solution. 
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INTRODUCTION 


The Jefimenko’s equations arise from an 
analogy between the laws that rules the 
electrodynamics and the laws of gravitation. 
Such analogy was proposed for the first time by 
Heaviside in a paper published [1] more than a 
century ago, where he supposed there must exist 
a second field due to moving masses and acting 
over moving masses only, called by Jefimenko, 
cogravitational field (sometimes this field is 
called Heaviside’s field). The Heaviside paper 
was forgotten for a long time until Jefimenko 
returned his work and made improvement to the 
Heaviside’s work in two books published and 
reissued since the 90’s decade [2], [3]. 


To start we need to write about the analogy 
made by Jefimenko between the laws of 
electromagnetism and the laws of gravitation. 
Oliver Heaviside proposed a system of 
equations with the same structure as the 
Maxwell’s equations, assuming the existence of 
a second field due and acting over moving 
masses only, called Heaviside’s field or co 
gravitational field and it is denoted by k, unlike 
the ordinary gravitational field g which is due 
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and acts not only on stationary masses but in 
movement. Although there are detractors of the 
Jefimenko’s theory of gravitation’ (see for 
example [5]), there exist books written by 
Wenceslao Segura [6] and Jolien E. D. Creighton 
[7] where they derived the Jefimenko’s equations 
from the linearized Einstein’s equations. Jefimenko 
derived the equations of gravitodynamics in a 
different way, starting from make an analogy of 
his retarded solutions of the electromagnetic 
field to the gravitodynamical field to finally get 
the analogous equations to the Maxwell’s 
equations for the gravitational and cogravitational 
fields. 


Next, we assume there is a configuration of the 
gravitodynamical field in an analogue way to 
those found in other works realized by 
Chubykalo and Espinoza [8], [9], where the 
authors have obtained the mathematical 
foundations on the Kapitsa’s hypothesis [10] 
about the origin of ball lightning related with 
interference processes. The configuration of the 
gravitodynamical field in our work begins with 


‘Called by us in a previous work gravitodynamical 
theory [4]. 
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the hypothesis that the gravitational waves exist. 
We are going to do a formal development about 
this issue based in the theoretical results 
obtained by Jefimenko. We will use the free 
gravitodynamical equations, that is, the 
equations valid for vacuum. 


It is important to emphasize that gravitational 
waves were predicted in the general theory of 
relativity by linearzing Einstein's field equations 
and this approximation is valid for weak fields. 
It is obvious that a complete theory of gravity is 
not linear but this linear approximation allows 
us to study a great variety of gravitational 
phenomena where gravitational induction is 
considered. In the same way that 
electromagnetic theory is not linear but Maxwell 
equations are applicable to a wide range of 
electromagnetic phenomena. 


Our results presented here do not try to replace 
the non-linear theories of gravitation, such as 
Einstein’s theory of general relativity [11] or 
Logunov’s relativistic theory of gravitation [12], 
but we want to show the importance of a linear 
gravitational theory, not only historical but 
methodological too, because we will show that 
exist properties of the weak fields that are not 
understood because they were not have the 
opportunity to appear in a linear theory of 
gravitation. 


THE  JEFIMENKO’S 
GRAVITATION 


EQUATIONS FOR 


The gravitodynamical theory is a generalization 
of the Newton’s gravitation theory, since the 
Newtonian theory of gravitation describes 
perfectly phenomena of a wide range of masses, 
but it does not consider the behavior of the 
fields of moving mass distributions. This is the 
reason why it is necessary to make an extension 
of the classical Newton’s theory of gravitation. 


We can write the Newton’s theory in a way as a 
force field theory in terms of the gravitational 
field g as 


V: g = —4ryg (1) 
and 
Vxg=0, (2) 


where y is the gravitational constant, ọ is the 
mass density given by ọ = dm/dV and dm is 
the element of mass contained in the volume 
element dV. 


As we know, the gravitational field g is defined 
by means of the force exerted by the mass 
distribution @ over a mass test 


2 


F, = mg. (3) 


In other words, the gravitational field is the 
perturbation of the space due to distribution of 
mass in some region which interact on a test 
mass m, Both masses (the mass creating the 
field and the test mass) can be moving or at rest. 


If we consider that the co gravitational field k 
exists, we need to define it in terms of the co 
gravitational force 


F, =m,(v xk), (4) 


Where F, is the force exerted by the 
cogravitational field over a test mass m, moving 
with velocity v. 


So, we can define the cogravitational field k as 
the perturbation of the space due to a moving 
mass distribution which interacts on a moving 
test mass. 


Jefimenko started to derive their gravitational 
equations from the next expressions 


enor | Pestle ser +f Por fer © 
and 
tE o 


where y is the gravitational constant, c is the 
velocity of propagation of the fields? and r is the 
vector directed from the element of volume dV’ 
(the source point) to the point where the 
gravitodynamical field is measured (the field 
point) and r is its magnitude, the square 
brackets designate that the quantities inside 
them are evaluated in the delayed timet’ = t — 
r/c. The integrals are evaluated over all space. 


We can see from (5) and (6) that the 
gravitodynamical fields have four causative 
sources, namely: the mass density ọ, the 
temporal derivative of the mass density 0,0, the 
mass current evand its time derivative 0; (eV). 


Jefimenko obtained the — gravitodynamical 
equations making use of the vector calculus and 
some vector identities. The gravitodynamical 
equations are: 


V: g = —4nyo, (7) 

V:-k=0, (8) 

T TE (9) 
ôt 


*Jefimenko assumed that the velocity of the 
propagation of the fields must bec, i.e, the finite 
speed of light. But we have demonstrated in [4] that 
this velocity must be finite or instantaneous. 
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where j= ọv is the mass current density. It is 
evident the analogy between the Maxwell 
equations and the Jefimenko’s ones is not 
perfect: while we have only one kind of mass, 
we have two types of electric charge. While the 
electric field is directed from the positive 
charges that generate it and is directed toward 
negative charges, the gravitational field is 
always directed to the masses by which it is 
created. Another difference is that the magnetic 
field is dextrorotatory (right hand) with respect 
to the electric current through which it is 
generated, while the gravitational field is always 
levorotatory (left hand) with respect to the mass 
current through which it is generated. In spite of 
these differences, the system of equations 
obtained by Jefimenko describe correctly the 
behavior of the weak gravitational fields, as we 
have already seen they are deduced from 
different formulations in [3] and in [6], [7]. 


We define a Gaussian system of units for the 
gravitodynamical field, in order to simplify our 
calculations. We will call this system 
Gravitational Gaussian System (GGS).To do 
this we need to introduce the next rationalized 
quantities: 


Tablel. Rationalized quantities in the new gravitational 
Gaussian system of units”. 


Gravitational field|Co gravitational field 
G=y 'g K=y'ck 
[G] = ML~* =Jef | [K] = ML~* =Jef 


Formula 
Units 


If we introduce such quantities in the system of 
Equations (7)-(10) we obtain the next system of 
equations in the GGS 


V-G=-—4n0, (11) 
V-K=0, (12) 
10K 
VxG= ar (13) 
10G AT, 
VXK= 3 ob (14) 


Where, j = ovis the mass current density and 
vis the velocity of the mass distribution 
generating the co gravitational field. 


GRAVITATIONAL WAVES 


The Jefimenko’s theory of gravitation predicts 
also the existence of gravitational waves. We 


*We define the unit Jefimenko abbreviated Jef for the 
rationalized gravitational and co gravitational fields 
in honor of Jefimenko. 
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Vxk-—=—=-— j, (10) 


will study the properties of such waves, in this 
and in other sections (especially in the section 
V, where we will study the energy and the 
Poynting vector of such waves). 


We will obtain in this section the wave equation 
for both fields, namely, the gravitational field 
and the cogravitational one, making direct 
calculations on the system of equations (11)- 
(14), we will see that these equations lead us to 
the wave equation. We start calculating the curl 
on the equation (13) 


10 
VxVxG=--—VxK, 15 
côt (5) 


and substituting in equation (14) and using the 
identity for a Laplacian of a vector, 
VxVxV=V(V:V)-—AV, (16) 


for any arbitrary vector V, and where A= Vĉ?is 
the Laplacian operator, we get 


10°G _ an (vo +42) z 

eoe ET ze) an 
the in homogeneous gravitational wave 
equation. 


In a similar way, starting from (14) taking the 
curl and using the identity (16), we get 
AK 10E ao, 18 
a T TVÍ (18) 
Is the inhomogeneous cogravitational wave 
equation. Both expressions (17) and (18) are 
field waves propagating on the space with a 
velocity c. 


If we consider regions without masses 
distributions and current masses we obtain the 
homogeneous wave equations, namely, 


10°G 
And 
1 0°K 


The equations obtained (19) and (20) can be 
solved by a sum of two vector functions, @, 
and Q2, 


Q,(k:r—ot)+@,(k:r+at), (21) 


Where k = (ky, ky, kz) is the wave vector’, 1 
and @ 2 are general expressions which represent 
plane waves propagating with velocity c in 


4 To avoid confusions, we use italic boldface k to 
represent the wave vector and normal boldface k to 
represent the cogravitational field. 
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opposite directions. Then the solution (21) for 
the wave equation can be derived from the 
method of separation of variables. If we 
introduce the harmonic dependence given by 


G = Geiko ng (22) 
and 
K = Kye! Tt ny, (23) 


Where ng is a unitary vector in the direction of 
propagation of G and nx is in the direction of 
propagation of K, the wave’s equation results in 
the dispersion’s relation 


An Interesting Wave Solution of the 
Jefimenko’s Equations for the Free Space 

The prediction of gravitational waves by means 
of the Jefimenko’s gravitodynamical theory let 
us search a various types of solutions for the 
gravitodynamical field in vacuum. For example, 
as we will show in this section, solutions exist 
that have a set of interesting properties. These 
solutions have the form of spheres of co 
gravitational field and ring-like gravitational 
field, in such way the total configuration of the 
gravitodynamical field oscillates. In such 
spheres the Heaviside’s field is tangent in all 
points over the surface of this sphere, and the 


In Electrodynamics is usual to refer to standard 
polarity in the solutions of the Maxwell 
equations when the electric field Eis a polar 
vector and the magnetic induction B is an axial 
vector, or pseudo-vector. This means that after a 
transformation of inversion of axis coordinates, 
E changes its signs, while B maintains its signs. 
Following the analogy between both theories we 
are going to look for solutions to these equations 
with standard polarity, to wit, when the vector 
G is polarand K is axial. We propose to solve the 
system of free Jefimenko’s equations by the 
method of separation of variables, we can write 
G and K as follows: 


GG, t) = y¥m)uO) (31) 
and 
K(r, t) = K(r)v(t), (32) 


where y(r )is a polar vector and K(r) is an 
axial one, also u(t) and v(t) are functions of 
time. 


Substituting (31) and (32) in the system (27)- 
(30) 


V-y(r) = 0, (33) 


w? — k?c? = 0. (24) 


For the election of sign k = + w/c, we obtain 


from the Jefimenko’s equations for 
gravitodynamical fields 
k-ng =Oand g-ny = 0, (25) 


@ 
(k x Ng) Go = — Konk and (k x Nx) Ko 
w 


From Eq s. (25) And (26) we can see that G and 
Kare vector mutually perpendicular to the 
direction of propagation and| Go| = |Kol. 


same for the ring-like configuration of the 
gravitational field. 


We will begin this section rewriting the 
Jefimenko’s equations for the gravitation, 
assuming there are regions of free space or 
vacuum, this means, ọ = 0 and j = 0. We get 


V:-G=0, (27) 
V:K=0, (28) 
VxG= ao (29) 
c Ot 
10G 
VxK=——.. (30) 
V-K(r) = 0, (34) 
1 1 dv(t) 
Vx y(r) = -zuo or “O (35) 
Vx k(r) = a GLY (36) 


We can equate the temporal parts of both 
equations to certain constants to obtain a 
consistent system, 

1 dv(t) a 
=Z es, 


And 


1 ut) 
vO at 
We are going to equating both equations to w, in 


order to obtain sinusoidal solutions and get only 
three constants in our system (37)-(38) 


(38) 


v(t) = Acos(wt — ô) (39) 
and 
u(t) = Asin(wt — ô), (40) 


where, A and 6 are arbitrary constants. 
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Thus, in this way, the equations for y and 
K become 


w 
Vxy(r) = zro), (41) 
and 

w 
Vxk(r)= 7 Y@). (42) 


Due to the linearity of the spatial components of 
the fields, we can add (41) and (42), then, we 
can define the vector 

o(r) = y(r) + x(r), (43) 


Such as, we obtain 
Vxo(r) = Zot). (44) 


First of all, we are going to solve (44), then, 
once we have obtained o, we can calculate 
y and K. We need to note that vector o has not 
polarity, but we can express its polar and axial 
parts as 


V(r) =5 fot) — o(-r)} (45) 


And 


r2 sin@ 00 0g rsin@ 


Gane 1 pea sing) _ seo E 1 pa 7 A(roy sin Na : spea E g 


We can obtain the next system of equations 
taking into account Eq. (44) and comparing it 
with (49) and (50) 


alog sin 0) _ woro, sin @ 


= - , (51) 
Orp) rog 

ðr c a 
and 
alra) Olo) _ ae (53) 


Or 00 c 


Expressing the variables 0,y og from (51) y (52) 
and replacing them in (53), we obtain the next 
partial differential equation for og, namely, 


1 272 


0 ; w 
and 90 (a4 sin 6)| + z Op 
=0. (54) 


a ð 
raaz% ) +o 


If we propose og = R(r)O(@) as a solution for 
(54) we find that R and © have to satisfy 


fe) + yr +A)}rR=0 55 
dr? c? er (2) 

And 

a 5 £ (Osi a} 10 =0 56 

dé \sin@ d8 a a (a6) 
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ap 


1 
k(r) = 5 {o(r) + o(=r)}. (46) 
Taking the curl of (45) and (46), and inverting 
the coordinates of equation (44), we have 
1 
Vxy(r) = z0 x o(r) — V x o(-r)} 


= Leo) + <o(-1)} 


A 
= zr) (47) 
And 
V x x(r) = Zv x o(r) + V x o(-r)} = 
12 olr) - 2o(=r)} = yn), (48) 


we can be sure that the system is satisfied. The 
only thing we need to do is to find the solution 
of (44), in order to find the solution of the 
system (41)-(42). To get such solution we will 
consider the vector o in spherical coordinates 
and we suppose that the solution has axial 
symmetry 


o = o, (r,0)° + og (r,0)ð + aop(r,0)@. (49) 
The curl of ø in spherical coordinates is 


(50) 


ðr r ðr a0 


where A is an arbitrary constant? If 2 = 0, then 
the solution for rR in Equation (55) must be 


wr » Wr 
A cos — + Bsin—, where A and B are 


constants, but, in general, A and B depend on r, 
so 


wr wr 
rR = A(r) cos —— + B(r) sin = (57) 


Now, we are going to substitute (57) in (55) and 
we obtain the next two equations 


ČA A, 2wdB a 
dr? r? c dr (38) 
And 

@B A, 2wda 5 
dr? r? cdr” (D 


Considering the fact that the coefficients of sine 
and cosine must be zero separately, due these 
functions has the same argument. 


To solve (58) and (59) we propose A(r) = ar™ 
andB(r) = br”, where the coefficients a and b 
are constants and n, m E€ N are constants, too. 
Then, we obtain the following characteristic 
equations, substituting the solutions proposed in 
Eqs. (58) and (59) 
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2w 
am(m — 1) + Aa + ma =0 (60) 
and 
2w 
bn(n — 1) + Ab — g anret =0. (61) 


Both equations can be satisfied for the next two 
cases: 
Casel) m = 0,n = —1,A = —2 and a = 
—wb/c, and regarding Eq. (57) for b = 1, we 
obtain 


1 ( wr wr 2 


R = — |- — cos— + sin—]; 62 
7z z cos sin— (62) 


Case2)m = —1,n = 0, A= —2 and b = —wa/ 
c, and regarding Eq. (57) for a = 1, we obtain, 


or . =). 


R=>( Z4 63 
=z cos z sinz (63) 


Then, from the solutions (62) and (63) we have 
the general solution for Eq. (55) 
C. 
R(r) = = (—= cos + sin) 
d i C y wr c wr . wr 
+ => (cos + —+—sin—), (64) 
Where Cı and C, are oe constants. This 


general solution can be expressed as 


and 


In order to write the solutions in a short way, we 
define 


wr wr | (wr 
Ç = cos (— - a) +—sin(—- a) 
c c c 


and 


z w*r? (es ) 
=— cos (— — q 
I c2 c : 


such as we can write the solution of Eq. (44) in 
spherical coordinates as 


o=é(- Z coso) +e (4 sind) 
+ë = sin 6) 0. (71) 
Where we have multiplied by é€w/c_ for 
convenience, andé has dimensions [č] = g cm. 


In this way, we can see that the component in 
@ direction corresponds to the vector y and the 
components in # and @ directions correspond to 
the vector K. That is, 


y= ¢(S sina) @ (72) 
and 


6 


R(r) = S {eos (= - a) +sin (= - a)}, (65) 
Where C and q are arbitrary constants. 


Now, we considering the Eq. (56), which for 
A = —2 become: 


d( 1 
70 alae Ade (Osin )} +20=0. (66) 


This equation has the general solution: 


O(6) = Cz sin@ + C,(coté 
— sin 8 In|csc8 —cot@|) (67) 


Where C3 and C, are arbitrary constants too. 


Due to in 0 = (2n + 1)m the corresponding 
solution has a singularity, we can make C, = 0. 
Also, due to the homogeneity of the equation for 
the vector o, we can make C3 = 1. 


In this way, we can write the solution of the Eq. 
(54) as follows 


04 (7,8) = 5 {cos (= — a) + = sin (= — a)} sin 8. (68) 


We can use the system (51)-(53) to find o, (r, 80) 
andøoọ (r, 0), namely 


o, (r,80) = sa {cos (= = a) + 2r sin (= = a)} cos@ = (69) 


Or 
cos (— - «)| sin 0. (70) 
k= e(Ž cosa) + ë (Zsino)a. (73) 


Finally, we write the solution for the 
gravitational and co gravitational fields as 
follows 


}sin(wt — ô) (74) 


a 
Il 
ww 
= 
Q 
£ 
3 N 
4. 
5 
D 
ns” 
© 


K= f (Bos o)? PHE (sin o) 8} cos(wt —6), (75) 
where we have bear in mind temporal solutions 
(39), (40) and the spatial solutions (72), (73). 


The necessary condition so that solutions (74) 
and (75) do not diverge in r = 0 is, 


wr wr | (wr 
£(0) = {cos (+ a)+ z sin (— a)} 
to fulfill such condition we need cos a = 0, this 


implies a = (n + 1/2)r, where n= 
0, +1, +2,... 


=0, 


r=0 


Now, we calculate the next limits to ensure that 
the solutions converge, 
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tie ie E E 
ge Wes = aes ea aes? 
these limits are evaluated for a = 1/2, also 
and 7 were expanded in power series of r. 


we have the energy density defined’ as 
w= ig +K?) (76) 
8r , 


so, we obtain the limits 
“3 
limG= 0; limK = coswtZ; limw 
r>0 r>0 r>0 


=a cos*wt, (77) 
Zis the unit vector in the direction Z + of the 
Cartesian system. 


We propose œ = mņt/2 and 6=0, because 6 
defines the initial wave phase of the fields G and 
K, we can write convergent solutions for these 
fields: 


G= f (5 ih 6) o} snöat (78) 
and 

K= f (= cos 6) PLE (3 sin 8) a} coswt, (79) 
where 

¢=—cos (“)+sin(*) and n =¢ - 


We can conclude that the solutions (78) and (79) 
for the Jefimenko’s gravitational equations for 
the free space involve the novel existence of 
spherical formations of the gravitodynamical 
field. 


Analysis of the Energy and the Energy Flow 
of the Gravitodynamical Field 


The expression of energy density given by Eq. 
(76) can be changed after some algebraic 
manipulations in another that contains a part 
time-dependent and an independent one, 
namely: 
& w2 2 
= meas 
2 2 2 
$ í uA +Tsine| 


167 || r6 rê 
272 


4 2 2 
sin? 0 + [5 cos? 0 + Tesina} 
r r 


“Tye sin? cos(2wt). (80) 


*See, page 303 of Gravitation and Co gravitation [3] 
by Jefimenko. But in our case we use GGS units. 
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From this expression we find the geometric 
places where the energy density do not depend 
on time t. Those geometric places are 


The points along the Z axis where is satisfied 


= , 


WZ WZ 
t pe 

an (—) - 
where 0 = 0, m and ¢ = 0. 


The surfaces where r satisfies 
wr? 

2 72 2 

n =ç ( zz — 4 cot o) 


The cross-section of such surfaces is drawn as 
discontinuous curves in the Fig. 3. 


Now, we are going to obtain the 
gravitodynamical energy Eg inside a sphere of 
radius R centered at the origin, by means of 
R n 2m 
Erg =Íf] w(r, 6, p, t)r? sin@ drd@dd = Eç(R) 
000 


+ E,(R,t), (81) 
where 
& [w*R* wR? wR 
Eg(R) = ana as 5 sin? ( i ) el, (82) 
and 
& 
Eg(R,t) = zp?" cos 2wt (83) 


We can see from Eq. (83) that gravitodynamical 
energy does not change in time within spheres 
of radiuses R which are solutions of the next 
equations obtained respectively making ¢ = 0 
and n = 0 


i (5 7 wR 94 
A aa) 
and 

wR ae 


The surfaces whose radii satisfy Eq. (84) 
contains only co gravitational field, there is not 
gravitational field over those surfaces, as we can 
verify from Eqs. (78) And (79) taking ¢ = 0. 


Now, we analyze the energy flow contained in 
the wave field given by (78) and (79). 


As a first step we are going to calculate the 
Poynting’s vector in GGSunits 
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S= KxG 
~ 4r 


_ & (won sin? @ 
~ 8r rm 

wg? sin(26) ~ 
E 2 Ja 


; | sin(2we), (86) 


Now, we will calculate the total momentum of 
the gravitodynamical field within a sphere of 
radius r centered at the origin. We will do this 
making use of the fact that the Poynting vector 
is proportional to the vector of the density of 
momentum, so that we can calculate the integral 
of the Poynting vector over the volume of the 
sphere of radius r. To make this easier we will 
express the unit vectors in spherical coordinates 
system in a Cartesian one, namely 


? = sind cos PX + sin 0 sing Ẹ + cos 0 Z And Ô = 
cos@ cos@X + cos @ sing y — sin 0 ĉ 


Integrating Eq. (81) over the given volume we 
obtain: 


I Sr? sin 8 drdðdo 


2w sin(2wt 2 
= a S sin‘@| dr2 
r 0 


= 0. (87) 


T 


We can interpret this result as follows. The total 
momentum of the gravitodynamical field (78)- 
(79) in a volume bounded by an arbitrary sphere 
centered at the origin is null at any time given. 


We will obtain the geometric places where the 
Poynting vector is zero at any instant. To do 
this, we need the conditions when the Poynting 
vector is zero and this is obtained by means of 
Eq. (76), 

Ç? sin20 = 0 and ¢nsin’@ = 0. (88) 
From the first equation of (88) we have 

¢ = 0, which satisfies both equations (88). We 
obtain the equation 

Wr Wr 
Accordingly, the geometric places for the case 
(1) are spheres whose radiuses satisfy Eq. (89). 


sin 29 = 0. Which means that @ can be 0, 2/2 
Orr. 


e 6=0,7, In this case both equations satisfy 
the conditions (88). Therefore, the 
geometric place is the Z axis. 


e @=n7/2, we have two possibilities to 
satisfy the conditions (88), namely: ¢ = 0, 


as in case (1) or 7 = 0. From this condition 


we have 
wr 
wr FA 
tan (—) = — (90) 
= 


So, we have the geometric places are rings in 
the plane z = 0 whose radiuses satisfying Eq. 
(90), these rings corresponding to the case 
0 =7/2 andņ = 0. In all points over these 
rings the cogravitational field is zero. 


The Poynting vector is tangential in all points 
over these surfaces. This can be seen from Eq. 
(81). This fact clarifies the conservation of 
energy within spheres of radiuses (85). 


The geometric places where the Poynting vector 
for the gravitodynamical field given by (78)- 
(79) is null at any time are: 


e Z axis, called cogravitational axis because 
the ordinary gravitational field does not 
exist there. 


e Rings at the plane z = 0 whose radiuses 
satisfy Eq. (90), called gravitational rings 
because there is not cogravitational field on 
them. 


e Spheres centered at the origin whose 
radiuses satisfy Eq. (89), called 
cogravitational spheres, because there is not 
gravitational field on them. 
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Figurel. Poynting vector field distribution of the 
gravitodynamical field at a given time in the plane 
2= 0, 


Let’s see the graph where is shown the 
distribution of the Poynting vector field in order 
to clarify the results obtained in this section. 
Due to the axial symmetry of the energy density 


°The graphics were performed in Mathematica™. 
Thez axis is the ordinate and the y axis is the 
abscissa, and we put in the program c = 1 and w = 1 
for simplicity. 


Open Access Journal of Physics V2e 11 e 2018 


Spherical and Ring-Like Configurations for the Gravitodynamical Field 


and the energy-flux density, we can consider 
only the distribution in the planex = 0. 


In Fig. 1 we can see the vertical cogravitational 
axis that matches with the z axis, And we ca 
also see the cross-section of three spheres, 
which we will call G-sphere to the first one, K- 
sphere to the second one, and G-sphere to the 
last one’, in an arbitrary instant. The total 
gravitodynamical energy conserves within G- 
spheres due to theenergy-flux vector at the 
surface of this sphere has tangential component 
only. 


We can see also that energy transfers from these 
G-spheres to the gravitational ring (the equator 
of such spheres) and after a period defined by 
the function sin(2wt) in Eq. (81) the movement 
is reversed. Inside the first G-sphere the energy 
transfers from the cogravitational axis to the 
gravitational ringand having spent some time 
returns. The energy within the K-sphere is also 
conserved, we can see this because the Poynting 
vector is zero in every point of the K-sphere 
graphed. The energy is transferred from the 
surface of the K-sphere to the gravitational rings 
of the G-spheres. An analogue exchange of 
energy occurs between next G-spheres and K- 
spheres. 


We want to emphasize the fact that the Poynting 
vector field reverses their direction after a time 
due to the function sin(2wt) present in Eq. (86). 


Let us see the cross-section of the Poynting 
vector field in the plane z = 0 in Fig. 2. 


Figure2. Pointing vector field distribution in the 
plane z = 0, for a given instant of time? 


And at last, we can see the graphic of the cross- 
section in the plane x = 0 of the surfaces where 
the energy density is constant, the graphic of the 


7 Called G-spheres because they have a gravitational 
ring in the equator, and K-spheres because they are 
spheres of co gravitational field. 

ŝXaxis is the ordinate and Y axis is the abscissa. 
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first K-sphere; the graphic of the second G- 
spheres. 
We need to note that these surfaces do not 


change in time in vacuum this means they do 
not deform nor displace when the time goes by. 
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Figure3. Cross-section of the surfaces where the 
energy density is constant (dashed lines). The 
continuous lines represent the K-sphere and the G- 
spheres. 


In Fig. 4 we obtain decreasing energy in the 
intervalR € [0,7], where we have plotted the 
graphics of the total gravitodynamical energy 
Erg for four different time values, namely, 
t= 0, 2/4, w/2 and 32/4. We have chosen 
these values due to the periodicity of the 
function cos 2wt in the time-dependent term of 
this energy. Here, we can see how the energy is 
changing in different the time values given 
before. There are various points where the 
gravitodynamical energy is constant for 
different values of time. For example, we can 
see that in the point (2.75, 3.8) all the curves 
intersect, this means, at R = 2.74 cm we obtain 
the total gravitodynamical energy Erg = 


—3.8 x 10-3erg?. 


Figure4. This shows how the total 


graph 
gravitodynamical energy Erç alternates as the 
distance varies for the different values of time time 
t = 0, 2/4, 2/2 and 3m/4 and for R € [0,7]. 


? These values of distance and energy are only for 
reference, because we have to remind that we have 
chosen the values w = 1 and c = 1. 
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Now, we want to show both dependences in 3D 
graphics, and we are going to analyze them. Due 
to the periodicity of the term time-dependent we 
will fix them for t € [0, 27]. First, we have Fig. 
5(a) the interval R € [0,1 ]. 


Figure5(a). Graph of the total gravitodynamical 
energy contained in the G- and K-spheres. Erg for 
the intervals R € [0,1] andt € [0, 27]. 


In the Fig. 5 (b) we show the total 
gravitodynamical energy in the intervalR € 
[0,10]. 


FigureS(b).Total gravitodynamical energy in the 
interval R € [0,10]. 
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Figure6.Contour plot of the total gravitational 
energy in the intervals R € [0,1.25]andt € [0, 27]. 
The cross-section where the total gravitodynamical 
energy is null, forms semi-ovoid. 


At last, we want to show the graphics of the 
cross-sections of both fields, in Fig. 7 we have 
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the gravitational field in the plane z = 0. In Fig. 
8 we have drawn cross-sections of the 
cogravitational field in the planes y= 0 and 
x = 0 respectively. 
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Figure7. Ring-like form of the gravitational field in 
the plane z = 0. 
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Figure8.Cogravitational field in the planesy =0 
and x = 0, this field does not have components in 
the plane z = 0. 


Convergent Solution (78) and _ (79) 


Represented as a Superposition of Two 
Divergent Solutions 


We can represent solutions (78) and (79) as a 
superposition of two waves spreading in 
opposite directions in each point, as the same 
way as Eq. (21). To do that only we need to do 
an algebraic transformation. 


We call 
Gravitational convergent solution. This G, is the 
superposition of the two waves G) and 


G()spreading in opposite directions at every 
point. In a similar way, we call 


co gravitational convergent solution. These 
solutions converge in r=0 & 6=(n+ 
1/2)x, where, n = 0, +1, +2, +3,... 
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Where in both cases we have 


čw sin@ wr wr | (wr pS 
G) = Zor? [cos( - wt) + z sin( F wt)| o, (93) 
čwsinð wr wr | (wr = 
Guy = Jcr? [cos ( P + wt) + sin = + wt)| o, (94) 
EcosO; or wr wr - 

Ko) = a [sin ( > wt) z cos( : wt) |? 

čsinð | or wr wr? wr ~ 
+ 73 | - os (= ot) +( = sin(—- ot) 0, (95) 

_ §cos@y wr wr wr 

Ki.) = ae [sin ( Z + wt) P cos( z + wt) |? 

€sin@[ wr wr wr*\ sar ~ 
+ 23 z cos( ; + wt) + (: -zzz }sin (—+ wt) 0. (96) 
The solutions (93)-(96) are solutions of the free Jefimenko’s equations. 
It is also possible demonstrate the next equations 

1 1 
and 
1 1 

K) = 2 (-K, + K.) and Kio) = 2 (K4 + K,), (98) 
where Gq and Kg are divergent solutions of the system (11)-(14): 
Gi = f (Ss sin o) 3 sin wt (99) 

£ cr? 
and 

26a a Na a 

Kı = f (Ecos o)? + § (sind) 8} cost, (100) 


where 


OA 
4 = cos (=) +“ sin(“) and Na = Cq — 2 sin( 


We said before that G5) ,G(. , Kí») and Ke) 
are solutions of the free Jefimenko’s equations 
for gravitation and they are divergent in r = 0. 
And we can conclude this section emphasize the 
fact that this kind of gravitational waves of 
Jefimenko’s solutions allows interference 
phenomenon as a superposition of two 


gravitational waves spreading in opposite 
directions. 
CONCLUSIONS 


As we have seen in section H, Eq. (21) 
represents a function of two waves spreading in 
opposite direction, both in k direction. Such 
waves are similar to electromagnetic waves, that 
is, G and K are transversal waves, perpendicular 
to the direction of propagation defined by the 
Pointing vector S. 


The superposition of the mentioned waves 
spreading one from the origin to infinity and the 
other one from the infinity to the origin 
produces stationary waves. Both waves have 
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T) 


axial symmetry. So, we have obtained a free 
stationary gravitodynamical field, it is 
consequence of gravitational interference 
processes. 


In this gravitodynamical configuration surfaces 
(the dashed lines in Fig. 3) and points (in the z 
axis) exist, where the energy density is constant. 
Such surfaces and points are nodes of energy 
density waves. 


The cogravitational spheres obtained, can be 
containers of mass, because if we consider Eq. 
(4), any cogravitational field exerts a perpendicular 
force on any mass to the plane formed by the 
vectors v and K, where v is the velocity of the 
particle and we have already seen that the 
cogravitational field in these unusual formations 
is tangent in every point of such surfaces. That 
is why the particles contained inside these 
spheres cannot leave such co gravitational 
spheres. The gravitational rings will take the 
particles of the gas contained in it and they will 
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turn them in a direction of rotation of such a 
field. 


From figure 4 we can see the points where the 
energy is constant, such points are those where 
the different curves are intersected and they 
represent the nodes of energy waves. 
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